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Abstract. - In the past twenty years, shear- banding flows have been probed by various tech- 
niques, such as rheometry, velocimetry and flow birefringence. In micellar solutions, many of the 
data collected exhibit unexplained spatio-temporal fluctuations. Recently, it has been suggested 
that those fluctuations originate from a purely elastic instability of the flow. In cylindrical Cou- 
ette geometry, the instability is reminiscent of the Taylor-like instability observed in viscoelastic 
polymer solutions. In this letter, we describe how the criterion for purely elastic Taylor-Couette 
instability should be adapted to shear-banding flows. We derive three categories of shear-banding 
flows with curved streamlines, depending on their stability. 
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"The stability of viscous liquids contained between two 
rotating cylinders" of radii Ri and i? Q -or Taylor-Couette 
(TC) flow-is the benchmark problem for instability of 
flows with curved streamlines. It was the title of a seminal 
paper by G.I. Taylor in 1923 [TJ, wherein Taylor showed 
that the purely annular flow eventually becomes unstable. 
Above a critical rotation speed, a secondary vortex flow 
sets in, with periodicity A *~ d along the vorticity direc- 
tion, where d = R a — Ri. The original study by Taylor 
concerned simple incompressible Newtonian fluids. But 
many fluids are non-Newtonian and exhibit viscoelastic 
contributions to the stress [2]. In 1990, Larson, Shaqfch 
and Muller showed that the TC flow of polymer solutions 
could also become unstable to a Taylor-like instability [3j . 
The kinematics of the unsteady flow are roughly similar to 
those of the Newtonian case, i.e. after a critical threshold, 
Taylor vortices appear, but the destabilizing mechanisms 
are very different, depending on two different kind of non- 
linearities. 

It is well known that Newtonian fluids can exhibit in- 
creasingly unstable flows for large values of the Reynolds 
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number. When only the inner cylinder is rotating, the 
Reynolds number depends on the rotation rate of the in- 
ner cylinder f^, such that Re = , where v is the 
kinematic viscosity of the fluid [II. In a simple Newto- 
nian fluid, the constitutive relation is the simple linear 
relation between stress and shear rate. Then, the only 
non-linearity in the equations of motion comes from the 
advective term on the velocity (v.V)v, in the equation of 
motion to ensure consistency between Eulerian and La- 
grangian descriptions of fluid motion. The Reynolds num- 
ber Re is linked to the relative magnitude of this term 
with respect to the dissipation terms [4]. 
In polymer solutions, and in many non-Newtonian fluids, 
the primary non-linearity usually comes from the consti- 
tutive relation rather than the momentum balance. The 
constitutive equation is dynamical, i.e. it relates to the 
stress relaxation dynamics and typically includes a con- 
verted derivative on the stress T 2 . In this convected 
derivative, consistency between Eulerian and Lagrangian 
descriptions requires a convective term, now applied on the 
stress (v.V)T, and material frame independence requires 
additional terms of similar dimensionality Vv.T 2j3 . The 
dimensionless group linked to the magnitude of those new 
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non-linear term is the Weissenberg number Wi = 7-7, 
where 7 = is the typical shear rate in the flow and 
r is the stress relaxation time [5] . The similarity between 
Wi and Re is more apparent by defining Re as a function 
of the viscous diffusion time t v( i = =y , Re = Tydj 1 6 . Re 
controls the magnitude of the 'inertial non-linearity , while 
Wi controls the magnitude of the 'elastic non-linearity' (7|. 
In general, both Re and Wi are non-zero, but in many 
practical cases for polymer solutions and melts, the elas- 
ticity number £ = ^ = is large, leading to negligible 
inertial effects. 

In the simplest TC flow, where only the inner cylinder is 
rotating, and in the small gap limit, i.e. d <C Ri, there 
exist two dimensionless groups, one relevant to the 'purely 
inertial TC instability' Sj = \/Ai?e-derived and observed 
by Taylor fTl-and one to the 'purely elastic TC instabil- 
ity' E e = yAW^-derived and observed by Larson et al. [3 . 



Here, A : 



Ri 



is the geometrical ratio linked to the stream- 



line curvature, necessarily finite for the instability to be 
linear [f][7]. Note that the Taylor number is usually de- 
fined as Ta = £ 2 |4j. In the purely inertial case, the flow 
becomes unstable for £i > to'. In the purely elastic case, 
the flow becomes unstable for S e > to. Both to' and to 
are coefficients of order unity, with precise values that de- 
pend on the boundary conditions [4j[8]. 
In this letter, we extend the expression of the in- 
stability criterion for viscoelastic 'shear-banding flows'. 
Shear-banding is yet another curious but ubiquitous phe- 



nomenon occurring in complex fluids 12|13 . When a fluid 



material is sheared, the strain rate gradient can be very 
large in narrow zones of the sample. Adjacent domains 
of markedly different strain rates are identifiable. This 
phenomenon has been observed in a variety of systems 
but in this letter, we focus on the steady shear-banding 
phenomenon of 'wormlike micelles' 12 13 . Entangled 
wormlike micellar solutions are model rheological fluids 
due to their linear Maxwellian behaviour for small defor- 
mations, characterized by a single relaxation time r and 
elastic modulus Go pi. Furthermore, the robustness of 
their shear-banding behaviour makes them attractive for 
the general study of banding phenomena in complex flu- 
ids 
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Roughly speaking, a shear-banding flow is remi- 
niscent of a first order phase transition. Above a lower crit- 
ical Weissenberg number Wi[, the shear stress plateaus. 
Then, until a second higher critical Weissenberg number 
Wih, the flow is inhomogeneous, split in two bands with 
local Weissenberg numbers Wii and Wih- To leading or- 
der, for Wi G [Wii,Wih], an increase in the value of the 
macroscopic Weissenberg Wi only increases the propor- 
tion a £ [0, 1] of the high Wi band, following a 'lever rule' 
Wi ~ aWih + (1 — a)Wi\ [12 13 . This scenario has been 



roughly confirmed experimentally with various techniques 
e.g. pure viscometry, velocimetry, birefringence, etc. 



13 



but many fluctuating behaviours were observed in all the 



gathered data 13 



In a series of recent experiments [15 -18 we have shown tices 



that the interface between the bands undulates due to an 
underlying secondary vortex flow that is mainly localised 
in the high Weissenberg number (Wih) band. Elastic in- 
stabilities similar to the one observed in polymer solu- 
tions could be the source of many of the observed spatio- 
temporal fluctuations. This rationale was reinforced by a 
recent linear stability analysis of the diffusive Johnson- 
Segalman (dJS) model 



19 , a viscoelastic 



constitutive 
In 
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model widely used to study shear-banding flows 
this letter we wish to rationalize experimental data on the 
shear-banding TC flow of wormlike micelles by determin- 
ing the appropriate form of the instability criterion in the 
case of shear-banding flow. 

Rheological and geometric scaling of purely elas- 
tic flow instabilities. — In the introduction, we dis- 
cussed the TC instability of the purely Newtonian fluid 
and the purely elastic fluid, which are both idealizations 
that facilitate our analysis, but which capture only the be- 
haviour of very specific fluids. In general, however, non- 
Newtonian fluids can exhibit other attributes such as a 
Newtonian solvent contribution to the stress, a spectrum 
of relaxation times instead of a single relaxation time r, 
and/or 'shear thinning', i.e. a decreasing viscosity with 
increasing shear rate |2j. Experiments conducted on such 
non-Newtonian fluids have documented the effects of such 
fluid rheology on the elastic TC instability [9]. To ratio- 
nalize these observations as well as to generalize the elastic 
instability criterion to different kinds of flows with curved 
streamlines, McKinley and coworkers established a gen- 
eral criterion for elastic instabilities [Tol[TT] . If Re ~ 0, 
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viscoelastic fluids are unstable if — > m , where 

x y ' 

= T xx — T yy is the first normal stress difference [2], 
is the shear stress, I is the characteristic distance over 



10 , and K 



which perturbations relax along a streamline 
is the characteristic radius of curvature of the streamlines. 
For a purely viscoelastic fluid, £ = Ut ~ FliRiT = Wi d, 
1Z ~ Ri and -ff 1 - — M 1 - ~ Wi 10 and we recover the 



for the purely elastic instability: 



Ni_ _ Ni 

T XH To 

criterion of Larson et al 
AWi 2 > to 2 <^ y/AWi > m 
In turn, framed with respect to the general criterion de- 
rived by McKinley et al. [l0|[TT], our goal is to deter- 

K 



mine the functional form of the dimensionless ratio t^ 3 - 



in terms of measurable quantities in the case of shear- 
banding flows. By analogy with polymer solutions, we 
would expect that this ratio can be expressed in terms 
of a relevant geometric ratio and an appropriately-defined 
Weissenberg number. 

Effective gap. — The relevant geometric ratio can in- 
deed be inferred from experiments through the notion of 
an effective gap. In our recent experiments 17p8 



ognized that the vortices were mainly localized in the high 
Wi band, and that each interfacial wavelength between 
the bands corresponded to a pair of counter-rotating vor- 
17 , as illustrated in Fig. [lj,. In our previous publica- 



tions, we had noticed that the wavelength increases upon 
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Fig. 1: Effective gap scaling, (a) Overlay of two visualization 
techniques showing the secondary vortex flow in the high Wi 
band for a ~ 0.4 [l7]. (b) Wavelength scaling, following A = 
nad, with n = 3.8 ±0.1. For a > 0.6, the spatio-temporal 
dynamics of the vortex flow do not allow us to extract a single 
wavelength [16|17 . For a < 0.05, the size of the band is smaller 
than our spatial resolution. 



increase of the global shear rate, so one could infer the scal- 
ing X/d ~ Wi 16 17 . Then, by combining this scaling 
and the lever rule we can establish that A = nad instead of 
A = nd, where n is a number of order unity, whose precise 



value depends on the boundary conditions. The extent of 
the high Wi band acts as the effective gap. Increasing the 
global Wi increases a and so increases A. The validity 
of this scaling is shown on Fig. |TJd by re-plotting 2X/d, 
i.e. twice the wavelength of vortices, against a instead of 
Wi ~ 



16 . 



Local Weissenberg number. — As explained in the 
introduction, in a shear-banding flow, the global value of 
Wi is not a good measure of the local Weissenberg num- 
ber in the parts of the flow that are unstable. Instead, the 
dimensionless group relevant to the flow instability is the 
local value of Wih in the high shear rate band. In the in- 
stability criterion, one must replace Wi by Wih- Accord- 
ingly, the criterion for elastic instabilities in shear-banding 
flows should involve the term 

S* = VaAWi h (1) 

It has been observed in experiments that increasing 
the concentration (c) of surfactant, or decreasing the 
temperature (8) tends to increase the value of Wih- 
This fact is illustrated in Fig. [2i in the flow curves of 



21.22 . Note that as 



two different surfactant systems 
increases, the dimensionless stress plateau decreases and 
its range of Weissenberg numbers increases. In particular, 
Wih shifts to higher values. For the most concentrated 
solutions, viscometric measurements had to be aborted 
because the sample was ejected from the rheometer. We 
believe that this phenomenon is due to an instability of 
the free surface of the system, driven by the underlying 
bulk viscoelastic instability. However, we also note that 
the instability of the free surface could be triggered by 
second normal stress differences 23 . From Eq. Q, 
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Fig. 2: Experimental and theoretical "flow-phase dia- 
grams" [2l]. (a) Open symbols: Measured dimensionless 
flow curves for varying [CTAB]=3, 7,10, 12, 15, 17,18, 22wt.% 
at fixed [NaN0a]=0.3M (replotted data adapted from Cap- 
pelaere et al. [22], permission from Springer). Closed 
symbols: Measured dimensionless flow curves for varying 
[CPCl+0.2NaSal]=2,4,6,8,10,12,21wt.% (courtesy of Berret et 
al. |21| ). The arrow points in the direction of higher c or lower 
9 [21] . In both cases, measurements were done using a cone- 
and-plate device. The flow curves of the two systems do not 
overlap, even when the stress and the shear-rate are scaled 
with Go and r respectively, i.e. in the framework of the dJS 
model, the two systems have a different value of the coefficient 
l a\ (b) Analytical dimensionless flow curves obtained for the 
dJS model in simple shear 
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The different flow curves are 
obtained for varying n. The color map gives the value of the 

scaled dimensionless criterion E^js = E<jj s \f^~i~ ■ The arrow 
points in the direction of lower n. 



note that solutions of high c or at low 6 are more likely 
to be unstable, owing to the larger values of Wih- 



The case of dJS. So far, we have suggested a 
new relevant dimensionless group for elastic instabilities 
in shear-banding flows, without appealing to any partic- 
ular rheological model. To reinforce our argument, we 
can investigate the form of the instability criterion for the 
diffusive Johnson-Segalman (dJS) model, which has been 
widely used to study shear-banding flows 



20 . Recently, 



it has even been used in numerical simulations confirming 
the presence of a secondary vortex flow triggered by a bulk 
viscoelastic instability in the high Wi band 19 . In this 



we 



model, the stress is taken as the sum of a 'polymeric' part 
T p and a 'solvent' part T s . The total viscosity of the fluid 
is the sum of a polymeric and solvent part t/q = rj p + rj s , 



P-3 



M.A. Fardin et al. 



with the zero shear rate value of the polymeric viscos- 
ity given by 77° = tGq. The polymeric stress varies non- 
monotonically with imposed shear rate and goes to zero at 
large Wi, such that rj a is the asymptotic value of viscosity 
for Wi — > 00. 

To evaluate the expression ^ , we need an analytic ex- 
pression for the stress ratio in the high Wi band. Let us 
symbolize this ratio by [■^ x -]h- In the small gap limit, we 
can assume that the stress profile across the gap is close 
to the profile in a plane Couette geometry. We can then 
use the inhomogeneous plane Couette solution recently de- 
rived by Sato et al. [24] . 

In plane Couette flow of the dJS model, it is common to 
express the total shear stress as T xy = 77=^=3 cr, and the 

first normal stress difference as Ni - 



VT- 
2G 



N 
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The 



parameter a is the 'slip parameter' of the dJS e quation . 
Shear-banding happens if \a\ ^ 1 and r\ = 3 



< 



20 
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N is a dimensionless normal stress difference and a is a di- 
mensionless total shear stress. In plane Couette flow, the 
momentum balance imposes that a is a constant across 
the gap, but N(y) is a function of the position y in the 
gap 24 . For steady flow in the shear-banding regime, 
Sato et al. established that 

\Ar 



r r 



N 



V2 
KS = K{a 



T)K) 



(2) 
(3) 



where K (y) = \/l — a 2 Wi(y) and S(y) = VI 
are respectively a dimensionless shear rate and a polymeric 
shear stress, both functions of the position in the gap. 
In dimensionless form, the addition of the polymeric and 
solvent shear stress is expressed by a = S(y)+r]K(y) [24] . 
In the shear-banding regime, Sato et al. have found an 
analytic solution for the dimensionless shear rate profile 
K(y) that follows a hyperbolic tangent profile between Ki 
and Kh [24] , with Ki < Kh given by 



Ki = 
K h = 



V2 



(4) 
(5) 



In the high shear rate band, K ~ Kh = VI — a 2 Wih- 
Thus from eqs. (2), ^ and ^ we can obtain the following 
expressions 



(6) 
(7) 
(8) 
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3 V VI — 2r?> 
Then, overall, if we set I ~ WihOtd and 1Z ~ Ri, we get 

E dJS = yfcRWi h f{ri) = E*/(»?) (9) 



Therefore, the result we obtain using dJS is slightly more 
complex than the naive criterion E* since it also de- 
pends on the viscosity ratio. For shear-banding we require 
■q < 1/8, so we have 0.7 < f{rj) < 1.3. This result is indeed 
not surprising, since we had obtained E* in analogy with 
the purely elastic case derived using the Upper Convected 
Maxwell model, where 77 = [3] . In the homogeneous 
and non shear-banding elastic case, adding a Newtonian 
solvent also modifies the dimensionless group by the addi- 
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tion of a function f*(rf) ~ y 
Note that the expression for E^jg can also be expressed in 
terms of the two dimensionless variables K and a. Indeed, 
from the lever rule, a = j^zzjt , an d from eq. (4) and (5), 
Ki and Kh can be expressed in terms of 77, which can be 
subsequently expressed in terms of a using eq. (2). Ulti- 
mately, one can reach the following equivalent alternative 
expression for E^js: 



A 



I -a 2 



MJS{K,<J) (10) 



where T, dJS (K, a) 



is a function 

of K and a only, whose precise functional form is a little 
too cumbersome to be written explicitly. Figure [2Jd. plots 
the flow curves computed from eqs. (2), (4) and (5) [24] , 
together with the magnitude of E^js- We can see that as 
the shear rate is increased, the proportion of the high Wi 
band increases, the magnitude of the scaled criterion E^js 
increases and the flow is increasingly prone to instability. 
By comparing the experimental flow curves in Fig. I|k and 
the flow curve derived in the case of dJS in Fig. [2}) one 
can see that the effect of decreasing the Newtonian solvent 
contribution 7/ to the total stress is similar to the effect of 
increasing the concentration of surfactant, or decreasing 
the temperature. 

Boundary conditions and classes of unstable 
shear-banding flows. — Generally, we expect the 
relevant dimensionless group for elastic instability in 
shear-banding flows to be E s f, = E*/*(?7), where /* 
is a function of the ratio between the zero shear and 
infinite shear viscosities. We expect the specific form of 
/* to depend on the constitutive model used to study 
shear-banding [12]. Elastic instabilities will generate 
a secondary vortex flow with wavelength A = nad for 
E s 5 > m. As mentioned already, the precise values of 
n and m depend on the boundary conditions. Of prime 
importance are the values of m obtained for 'soft' (m s ) 
or for 'hard' (rrih) boundary conditions 18 . Essentially, 
the 'hard' case usually corresponds to a no-slip Dirichlet 
boundary condition, while the 'soft' case usually corre- 
sponds to imposing only continuity of the stress, i.e. a 
Neumann boundary condition. In both the purely inertial 
case [2] and the purely elastic case [8], it is known that 
m s < rnh- For a banded flow with Wi € [Wii, Wih), the 
interface with the low Wi band acts as a soft boundary 
for the high Wi band. But for Wi ^ Wih, a — 1, the flow 



p-4 



Criterion for purely elastic Taylor-Couette instability in the flows of shear-banding fluids 




0.001 



100 



Max[Wi h ,Wi] 



Fig. 3: Schematic instability diagram in the plane 
(Max[Wih,Wi],aA). The black lines represent the stability 
limits for soft and hard boundaries, £ s (, = m 3 <S> qA = 



(m s /Wih) and E s 



<4> aA = (m h /Wih) 2 , where we 



have arbitrarily chosen m s — 1 and nih = 3. The dashed black 
line represent the value of 1 x A = 1.13/13.33, the maximum 
curvature corresponding to our recent experiments [15}{18] . 
Above this line, the shaded region is inaccessible. The three 
paths 1., 2. and 3. illustrate the three possible types of shear- 
banding. The direction of the arrows represent the path fol- 
lowed by the state of the flow as the global Weissenberg Wi is 
increased. a c , a c i and a C 2 are the critical proportions of the 
high Wi band at which the flow state crosses a stability limit. 
Wi c is the threshold at which the type 1. trajectory becomes 
unstable for the first time, and at which the type 2. trajectory 
becomes unstable after a short relaminarization. 



becomes homogeneous again and the boundary switches 
from soft to hard. 

Therefore, for a given geometry, i.e. a given value of A, 
we can use basic Boolean logic to classify shear-banding 
flows into three possible categories depending only on the 
value of Wih- 

1. For sufficiently low Wih~i.e. high 9 and low c-the 
shear-banding flow is stable for any a, since E s b < m s 
even for a = 1. The flow can then become unstable for 
Weissenberg numbers above a critical value Wi c > Wih 
as in the case of a regular viscoelastic fluid, i.e. following 
the scaling Ej = \/AWi. 

2. For intermediate values of Wih-i.e. intermediate 9 and 
c-the shear-banding flow is unstable above a critical value 
a c when E s t> > m s for a > a c . Then as the imposed shear 
rate is increased and a — > 1 the boundary conditions 
change and the flow is stabilized, because the flow is below 
the threshold m/,. Eventually for Wi > Wi c > Wih the 
flow becomes unstable again. This case was the one we 



characteristics may change 

The three possible shear-banding scenarios can be 
illustrated on a stability diagram in the plane 
(Ma,x[Wih,Wi],aA), as presented in Fig. [3j When 
the global Weissenberg number Wi is increased above 
Wii, the flow state is given by a constant abscissa 
depending on the value of Wih (which is a function of 
the concentration and temperature of the solution). As 
Wi increases, the thickness of the high shear rate band a 
increases and so the state of the flow moves vertically to 
larger ordinates. Once the entire gap is filled, a A reaches 
its maximum depending on the geometry of the chosen 
TC system. Then, since Wi > Wih, the state of the flow 
is given by a constant ordinate A and moves horizontally 
as Wi increases. Any flow state with aA < A will be 
stable if below the stability limit E s & = m s , and unstable 
if above E s f, = m s and a fortiori if above E s t = mh- Any 
flow state with aA = A will be stable if on the left of the 
stability limit E s {, = wih, and unstable otherwise. 

Interaction with interface modes. So far, we 
have only considered elastic instabilities arising in the bulk 
of the high Wi band. But there exist other elastic insta- 
bility mechanisms 
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In particular, Fielding has shown 
that the jump in normal stresses between the bands could 
generate interfacial modes, even in plane Couette flow 20 



In her recent study in TC flow, Fielding suggested that the 
interfacial and bulk elastic modes lie in two separate re- 
gions of the space (A,7Vi|/j), i.e. of the space (A,Wih) 19 . 
The bulk mode prevails at high Wih and high curvature 
A. The interfacial mode prevails at low Wih and low A. 
Fielding's study would suggest the existence of another 
unstable region in the lower left corner of the stability di- 
agram sketched in Fig. [3j Nonetheless, only axisymmetric 
perturbations were considered in Fielding's study 19 , and 



the stability analysis was performed for a single value of a 
and rj. Interfacial and bulk modes may actually interact 
through non-axisymmetric mechanisms [26] , 

Wall slip and non-local effects. — We believe that 
the instability criterion we have derived for shear-banding 
flows can be a powerful guide to interpret experiments on 
wormlike micelles. Nonetheless, the criterion is fallible. 
In particular, we think that two additional phenomena 
can strongly compromise the validity of our scaling, 
since both have been shown to be relevant in some 
experimental situations. In both phenomena the local 
Weissenberg value in the high shear-rate band may not 
be equal to the upper boundary of the shear-banding 
regime on the flow curve. The first phenomenon is wall 
slip, which has been reported recently and may actually 
be a common feature of many shear-banding flows 
The second phenomenon is geometric confinement. 
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present scaling may be inadequate if 'non-local effects' 
become dominant 



28 



observed in our recent experiments 18 . 
3. Finally, if Wih is high enough-z.e. for low 9 and 
high c-we have two critical band widths a c \ and a C 2- 
For a > a c \, E s fc > m s . And for a > a C 2, S s b > wih- 
In this case, there is no stabilization for Wi > Wih- 

The flow remains unstable, although the spatiotemporal stress diffusion coefficient 20 24 



The 



'Non-local effects' are apparent in 
confined geometries when the size d becomes comparable 
to the typical interfacial width £ ~ fim, linked to the 

Even in a macroscopic 
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geometry with d 3> £, non-local effects can be important 
when the lateral extent of one of the bands is very small, 
i.e. a ~ or a ~ 1. Those effects were ignored in the 
analytic solution for dJS proposed by Sato et al. but can 
actually be derived directly from the dJS equations |29| . 

In summary, we have derived a useful dimensionless cri- 
terion to rationalize the onset of secondary flows in the 
base shear-banding flow of wormlike micelles. The valid- 
ity of the criterion for the case of dJS could be checked by 
numerical simulations for various value of the solvent ratio 
r\, and for a range of gap spacings (A) and Weissenberg 
numbers. On the experimental side, we are currently un- 
dertaking a large study of the stability of shear-banding 
flows for many different surfactant types, concentrations 
and temperatures. Preliminary results confirm the exis- 
tence of the three distinct scenarios that we derived here. 
Ultimately, the criterion could be extended to other flows 
with curved streamlines, if the localization and number of 
bands is known. 
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